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4 : Abstract 

Q . The influence of small additive noise on structure formation near a forwards and near an inverted 

bifurcation as described by a cubic and quintic Ginzburg Landau amplitude equation, respectively, 
O 1' is studied numerically for group velocities in the vicinity of the convective-absolute instability 

where the deterministic front dynamics would empty the system. 
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I. INTRODUCTION 

The formation of macroscopic structures in systems that are driven out of thermal 
equihbrium by an externally imposed generalized stress are usually investigated by deter- 
ministic field equations. However, under specific circumstances the influence of external 
deterministic or stochastic perturbations and of internal thermal noise on the pattern for- 
mation process should be taken into account to achieve a more realistic and quantitative 
description of experiments. One prominent example are the so-called noise sustained struc- 



turesl2, 



^, ^, ^ „^ in the convectively unstable parameter regime 

Q i, e.., t»-Coue„e ffflfl B Q, the RayM.h-B.na.d fl fl Q =..e.. 
or nonlinear optics [1^. Further examples are certain open-flow instabilities , e.g., in wakes 
and jets that are reviewed in jioj l. 

The noise sustained structures arise when an 

externally imposed through-flow or an internally generated group velocity v is large enough 
to "blow" the pattern out of the system according to the deterministic field equations. In 
this driving regime one observes in experiments 0,0,1^, 12, 1^ structures that are sustained 
by sources that generate perturbations in the band of modes that are amplified according 
to the supercritical deterministic growth dynamics in downstream direction sufficiently far 
away from the inlet. 

The criterion 3, at which f the pattern is blown out of the system under determin- 
istic laws which gave the threshold for the appearance of the noise sustained, supercritically 
bifurcating patterns in the above described experiments is a linear one. It was nonlinearly 
extended by Chomaz to the question of the propagation direction of nonlinear determin- 
istic fronts in infinite systems that connect the unstructured state to the finite-amplitude 
structured one. 

Here we study and compare the noise sensitivity of pattern forming systems in which the 
above described fronts are linear or nonlinear ones. To that end we investigate the cubic 
Ginzburg-Landau amplitude equation (GLE) for a supercritical forwards bifurcation and the 
quintic GLE for a subcritical inverted bifurcation, respectively, in one spatial dimension. 

We solve the GLE with additive stochastic forcing numerically. Our systems are finite but 
sufficiently long to allow the establishment of a statistically stationary large-amplitude bulk 
part - provided the latter is possible with the boundary condition of vanishing amplitude 
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at the ends. We focus our attention to parameters in the vicinity of the convective-absolute 
threshold at which the fronts of the deterministic GLE cease to propagate. And we investi- 
gate in particular the statistical dynamics of phase and amplitude fluctuations in the front 
region. 

II. SYSTEM 

We consider the stochastic, ID Ginzburg-Landau equation 

(dt + vdM = if^ + dl+ gMl'' + g,\A\^)A + ar] (2.1) 
for the complex amplitude 

A = ^A + i'^A = Re'^ (2.2) 

depending on x, t. Here 3fJ(Q') denotes the real (imaginary) part and R = \A\ is the modulus 
and $ the phase of A. The coefficients in (12.111 are taken as real for simplicity. We checked 
however that taking into account the (small) imaginary parts, that appear e.g. in the case of 
transverse Rayleigh-Benard convection rolls propagating downstream in a small externally 
imposed lateral through-flow or in the case of downstream propagating Taylor vortices 
P-llS^ does not change the major findings presented in this paper significantly. We consider 
the group- or mean fiow velocity f > in positive x-direction and the linear growth rate n 
of A as control parameters. 

We investigate two fixed combinations of the nonlinear coefficients {gziQb) that we refer 
to in this paper as follows 

(^3 = — 1 , (^5 = : cubic GLE (2.3a) 
(73 = 1 , (yf5 = — 1 : quintic GLE . (2.3b) 

The quantity a in measures the real strength of the complex stochastic force 

r/(x,t) = 3?r/(x,t) +i3r7(a;,t) (2.4) 

with statistically independent real and imaginary parts Sftr/ and '^rj, respectively. Both are 
Gaussian distributed with zero mean and 5-correlated such that 



< r]{xiM)Hx2M)]* > = 25{xi- X2)5{ti-t2) 



(2.5) 
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A. Unforced homogeneous solution 

We are interested in the effect of small additive noise on the spatio-temporal structure 
formation in large but finite or semiinfinite systems. Nevertheless it is useful to briefly recall 
first the properties of the most simple solutions of the unforced GLE in an infinite system. 
This shows what one might expect to see in the bulk of a very large system far away from 
the boundaries — ignoring for the moment questions related to boundary induced pattern 
selection processes. 

The GLE ()2.1|) shows for cr=0 a continuous family of traveling wave (TW) solutions 

A{x,t) = Re''^^''-'"^ (2.6) 

with constant wave number q, frequency Q = —qv, and modulus R given by 

^i-q^ + gsR^ + g,R^ = 0. (2.7) 

This TW solution family bifurcates at the marginal stability curve, fi = q'^, of the A=0 
solution out of the latter while the former becomes unstable there. The critical values are 
f^c = Qc = flc=0. The bifurcation is nonhysteretic and forwards in the cubic case 

i?2 = - g2 (2.8) 

and hysteretic, backwards in the quintic case 

R' = \±^fi-q' + \. (2.9) 

Here the lower sign refers to the lower unstable TW solution branch that exists for — | < 
/i — < 0. The upper TW solution branch identified by the + sign in Eq. ()2.9|) exists 
beyond the saddle- node bifurcation value n = q^ — \- These TW solutions are stable for 
wave numbers outside the Eckhaus unstable band [l9|. 



B. Convective-absolute instability 



The noise susceptibi 



ity of the pattern formation process described by the GLE ()2.ip 



changes significantly % |l6[ when crossing the parameter combination of /x, v shown in Fig. ^ 
for the so called convective-absolute instability This combination 

cubic GLE 



_3_ 

16 



v"^ + -^v - 1 ) quintic GLE 



(2.10) 
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is marked by the front solution of the deterministic GLE with a = undergoing a reversal 
of the front propagating direction in an infinite system. Consider a front that connects 
the basic state A = being realized at a; ^ — oo to a homogeneous solution with A 
at X — s> oo. For parameter values below (above) the respective curves in Fig. Qthis front 
moves to the right (left). Thus the basic state A = (the homogeneous solution A ^ 0) 
expands to the right (left). The region below (above) the respective curves in Fig. ^ where 
the basic state A = (the homogeneous state A ^ 0) invades the whole system is called the 

nn 

convectively (absolutely) unstable region of the A = solution [2, UM ■ Thus, the boundary 
fl2.1()|l is also called the convective-absolute instability boundary. ^ 

For the cubic GLE the boundary fic-a = results from a linear analysis 2]. For the 
backwards bifurcating solution arising in t he q uintic GLE the respective front that reverts 
its propagation direction is a nonlinear one |2C| . Note that in the latter case the convective- 



absolute instability boundary [17| connects for f ^ to the so-called Maxwell point ^c-a 
jiM = —3/16: For this value the minima of the potential V{A) = —^A^ — + ^A^ have 
equal height V = 

The boundary condition A{x = 0,t) = that we apply in our simulations stops any 
front propagating to the left and it changes, i.e., it deforms the front profile when the 
front is sufficiently close to the boundary at x = 0. This can be seen in Fig. |21 for the 
example of the deterministic quintic GLE. There the lines show the modulus profile R and 
the spatial growth rate k = R' / R versus x together with k versus R obtained numerically for 
several parameter values above the convective-absolute instability boundary. To facilitate 
comparison of different cases we introduce the reduced horizontal "distance" 

A = ^-1 (2.11) 

Vc-a 

from the boundaries shown in Fig. ^ Here 

f2,/u cubic GLE 



1 + 4/i) ~ Y 3 quintic GLE 

denotes the convective-absolute instability boundary ()2.10p . 

The results that we present here were obtained for > 0, i.e., in a situation where 
the basic state A = is unstable. For the backwards bifurcation in the quintic GLE with 
negative growth rates — ^ < /i < for which the above cited potential has a minimum 
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ai A = the situation is more complicated Q]: Not only does the establishment of the 
final front connecting the inlet condition A = with a statistically stationary saturated 
bulk with \A\ = 0{1) depend sensitively on the initial condition [say, A{x,t = 0) = 
versus \A\ = 0{1)] in the absolutely unstable regime, A < 0. But more importantly, in the 
convectively unstable regime, A > 0, we found that small noise does not seem to be able 
to generate with the boundary condition A{x = 0) = a noise sustained finite-amplitude 
structure with (|v4p) of order one when fi < 0: The deterministic front dynamics drives 
the large-amplitude part downstream and eventually any finite system is filled only with 
small-amplitude fiuctuations of A around the stable fixed point A = of the unforced 
system. 



C. Noise strength 

For the quintic GLE we choose the noise strength a = 10~^. The noise intensity should 
be compared with the minimum of the potential 

V{A) = -f^A' - gs^A^ - g,^-A' . (2.13) 

For our quintic case {g^ = 1, g^ = —1) the minimum at = Rj^ = | + is V{Rn) = 

-2^[l + 6/i + (l + 4/i)3/2]. Thus, the noise "temperature" cr^ measured in units of V{Rn) is 
a'^ /\V{Rn)\ = 9.2 10~^ for the control parameter = 0.05 that we have used in most of our 
calculations. 

A rough estimate for an equivalent noise strength for the cubic GLE would be to demand 
that the reduced noise "temperature" a'^/V{Rj\f) is in both cases the same. This would 
require for the cubic GLE at a common /i of, say, 0.05 that a is by about a factor of 13 
smaller than for the quintic GLE. 

However, basing the comparison on the requirement that a'^/V{RN) is the same for the 
cubic and quintic case one has to keep in mind that the curvatures of V around the states 
A = and A = Rn which are connected by the fronts remain different - cf. Fig. El Since 
these curvatures around A = {A = Rjy) measure the growth (decay) rates of fiuctuations 
around the respective states it is useful to compare their ratios via a kind of Ginzburg 
number G = \ V" {0)\/V" (Rn)- One has G cubic = 1/2 independent of /z and gs while Gquintic = 
/i [l + 4/i + y/1 + 4/i] ^ with /2 = —g^fJ'/g'i- Thus for fi = 0.05 and g^ = l,g^ = —1 one has 
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Gcubic — 23 Gquintic- This largely explains the stronger noise sensitivity of the cubic GLE for 
our parameters. In view of it we investigated the whole range of a between 10~^ and 10~^ 
for the cubic GLE. 

The cubic GLE with additional (but very small) complex coefficients has previously been 
investigated, e.g., for noise strengths of about cr = 1.9 ■ 10~^ in our units of eqs. ()2.HI2.5j] . 
The corresponding noise "temperature" a^/V{R]\f) is about 10~^ for a typical value of, say, 
/i = 0.035 This noise was found to fit the experimental results on the noise sustained 
traveling Taylor vortices under statistically stationary fronts in the convectively unstable 
regime of open Taylor-Couette systems with axial through-flow j^. 



D. Numerical methods 

Equation 1)2.11) was solved numerically with a forward-time, centered-space method |2ll 
subject to the boundary conditions 

A{x = 0,t) = = A{x = L,t) (2.14) 

on the complex amplitude. System sizes L were chosen to be sufficiently large to allow for 
the establishment of a saturated bulk amplitude. Typically, a spatial step dx = 0.4 was used 
with a time step of dt = 0.072. Calculations were performed for sequences of the paramater v 
at several values of the control parameter fi. Most of them were done at = 0.05. The noise 
source rj was realized by Gaussian distributed random numbers of unit variance that were 



divided by V dtdx to ensure independence of the correlation functions of the discretization. 
A test of different pseudo random number generators, namely, L'Ecuyer's method with Bays- 
Durham shuffie |21,], ranS |2l|, and the R250 shift-register random number generator 



22 



gave similar results. 

After the simulations were started, a sufficiently long time depending on the parameters, 
e.g., on the closeness to the convective-absolute threshold had to be waited until the system 
relaxed into a statistically stationary state with time independent averages. Thereafter 
time averages were evaluated over several consecutive time intervals and finally averaged. 
Within the forward-time integration method A{x,t) remains uncorrelated with r]{x',t) at 
the same time, < f{ri)g{A) > = < firf) >< g{A) >, so that, e.g., < Arj > = as well as 
< qr] > = 0. But < Vtr] > 7^ 0. Here the frequency Vt (wave number q) is defined as a 
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forward-time (centered-space) difference of the phase ()3.2|) . 
III. RESULTS 

The influence of additive noise on the pattern formation process described by the GLE 
()2.1|) is described in this section. 

A. Growth length i 

In Fig. m we show how the growth length i of the downstream pattern occurring in the 
forced cubic GLE varies with noise strength a. Here i is defined by the distance from a; = 
at which the root-mean square y/^^Ap^ of the fluctuating complex amplitude A reaches 
half its bulk value. In the absence of noise i diverges at the convective-absolute threshold 
V = Vc-a since there the deterministic pattern is blown out of the system. 

For finite a the solution with finite A is noise- sustained in the convectively unstable regime 
A > j^. In this regime i is far from the convective-absolute threshold well described by 
the relation i ~ —(1 -|- V2A) In a following from a quasilinear analysis of the cubic GLE 
presented here in an appendix. However in the vicinity of the threshold A = the growth 
length i obtained from the nonlinear GLE shows a characteristic crossover to the behavior 
at A < 0. 

The noise influences also in this absolutely unstable regime, A < 0, the finite amplitude 
solution at least close to threshold: The curves i{A, a) in Fig. EJ^a) break away from the 
dotted £{A, cr = 0) reference growth length curve at negative A values that decrease with 
increasing a, i.e., further and further away from the convective-absolute threshold. The 
associated inflection points can be most easily identified by the maxima in di{A,a)/dA 
shown in Fig. Efb). These peak positions of di/dA vary with a as shown in the inset of 
Fig. El^b) . So the growth length shows for the cubic GLE a definite noise sensitivity also in 
the absolutely unstable regime. 

This sensitivity is significantly smaller in the quintic GLE. This can be seen by comparing 
the behavior of the growth length with the fluctuations of the modulus R = \A\, of the 
frequency, and of the wave-number (cf. Sec. IIIIB|) . To that end we show in Figs. Eland 
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and di/dA together with the inverse of the standard deviations of the modulus 



sr = V< > - < R>\ 



(3.1) 



of the frequency sq (j3.9p . and of the wave-number Sg fj3.9|) at ;U=0.05 as functions of A 
for the cubic and quintic GLE, respectively. The noise strengths a = 2.5 ■ 10~^ and 10~^, 
respectively, used for these figures are roughly equivalent based on the criterion described in 
Sec. Ill (]l However, the potential minima in the cubic case are broader than in the quintic 
case - cf. Fig. 01- and therefore the modulus fluctuations in the former are larger than those 
in the latter one. This can be seen by comparing the reduced inverse y<B?> / Sr in the 
absolutely unstable regime, A < 0, of Figs. Efc) andEfc). 

The peak position of dl/d/S. coincides with the drop-off in the inverse standard deviations 
1/s. For the cubic GLE (Fig. it occurs at A=-0.049, thus being shifted significantly into 
the absolutely unstable regime while that of the quintic GLE (Fig. ^ remains at A=0. 

As an aside we mention that for the quintic GLE at a subcritical growth parameter of, 
say, /i=-0.05 the behavior of the growth length £ and of di/dA is for A < similar to the 
one shown in Fig. IHl^d) for /i=0.05. For /i < 0, A > we did not find a noise sustained 
large-amplitude solution. 

B. Frequency and wave-number correlations 

Previous investigations of the forced cubic GLE in the bulk part of the solution at far 
downstream locations x ^ i showed for different but small noise strengths that frequency 
fluctuations are in the absolutely unstable regime much smaller than in the convectively 
unstable regime i3|. In order to study this question of the noise sensitivity in both regimes we 
have investigated in more detail the frequency and wave- number fluctuations ai x = i/2,i, 
and 2i. The results are shown in Fig.Elfor the cubic GLE and in Fig.lHlfor the quintic GLE. 
Before we discuss them we first present some basic properties of the phase fluctuations as 
described by the forced GLE 1)2.11] . 

The phase fluctuations $ of the complex amplitude ()2.2j) define the frequency Q and the 
wave number q 
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respectively. Here dot (prime) denotes temporal (spatial) derivative. The growth rate k, of 
the modulus is given by 

= ^ = ^ V (3-3) 



R 

By means of Eq. ()2.1|1 the frequency can be expressed as 

n = i2K - v)q + q' + {7]A*) . (3.4) 

R"^ 

This relation holds for the cubic as well as for the quintic GLE with real coefficients. By 
squaring and averaging Eq. ()3.4|1 one gets the correlation functions 

< > + < > +2v < fig > + < g' ^ > -2 < fig' > -2v < qq' > 
- 4v < nq^ > +4 < K^g^ > -4 < /tfig > +4 < /tgg > 

(3.5) 



2 < > 



On the r.h.s. we have used the fact that within our forward-time integration method A 
remains uncorrelated with rj at the same time and we have approximated < 1/R^ > by 
1/ < i?2 >. 

Given that < \ri(t,x)\'^ >= 2/ dxdt in our finite difference simulation it is convenient to 
scale all correlations in Eq. (j3.5p by the quantity 

cr^ 1 o cubic GLE 

S2 = , Rl={ I (3.6) 

K dxdt i + + 1 quintic GLE 



thereby removing the singularities from the reduced correlation functions. For example one 
finds that 

<{n + vq-qr> ^ R% 7^ 

Here we have neglected the second line in Eq. ()3.5j) since all correlations in Eq. ()3.5|1 involving 
the growth rate k, are very small. 

< fi^ > is typically two orders of magnitude larger than < g^ > in the absolutely unstable 
regime, A < 0, - cf. Figs. El and IHl discussed further below. There the only contributions 
to Eqs. ()3.5I3.7|) of the same order as < fi^ > are < fig' > and < g'^ > - all the other 
correlations can be neglected - and furthermore < fig' >~< g'^ >. Thus, 

< fi2 > ~ < g'2 > (3.8) 
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in the bulk part of the system with saturated amphtude where < i?^ >~ Rj^. However, in 
the convectively unstable regime, A > 0, with much larger phase fluctuations the situation 
is more complex. Here < > is larger than < > except for the upstream region where 
the reverse holds. 

In Fig. and Fig. El we show the inverse of the standard deviations 



V<n'^ > - <n>^ , Sq = V< g2 > _ < g >2 ^ (3.9) 



reduced by S flH.(j|l for the cubic and quintic GLE, respectively, as functions of A for a; = 
i/2,i, and 2i. For the parameters shown in Fig. [SJand Fig. IHlthe mean frequency < Q > 
as well as the mean wave number < q > are negligible. Plotting the inverse of sq, Sg, 
and sr allows to visualize the small fluctuations in the absolutely unstable regime better 
than in a direct p.ot of. say. 4- SuA p.ots for 4 have been presented previously for the 
small noise strengths occurring in Taylor-Couette experiments [3]. On the lower level of 
resolution inherent in this data presentation these results show similar behavior as ours. 
However, plotting 1/sq instead allows to identify more clearly the crossover behavior from 
the parameter regime with small fluctuations to the one with large ones. 

The A- variations of I/sq, 1/sq, I/sr, and of di/dA indicate that this transition is shifted 
to negative A, i.e. into the absolutely unstable regime. A similar result for the transition 
between deterministic and noise sustained standing wave solutions of complex coupled cubic 
GLE's was deduced from the behavior of the second moments of the frequency and wave- 
number power spectra of the fluctuating amplitudes : With decreasing /i the correlation 
length defined via the time average of the second moment of the Fourier spectrum of A{k, t) 
begins to decrease towards values characteristic for noise-sustained structures in the con- 
vectively unstable regime clearly before /ic-a is reached when noise is present. Similarly the 
width of the frequency power spectrum starts to increase with decreasing fi already above 
the convective-absolute threshold fic-a [10]. 

However, the variation of 1/sq with A shows for the cubic case in Fig. El a broader 
crossover interval between large frequency fluctuations in the convectively unstable regime 
at A > and small frequency fluctuations in the absolutely unstable regime at A < than 
the curves and 1/sr for wave- number and modulus fluctuations. The A- value at which 
1/sg and 1/sr drop down towards zero agrees quite well with the peak location of di/dA. 
The latter moves with increasing noise strength further into the absolutely unstable regime 
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as shown, e.g., for the cubic GLE in the inset of Fig. Efb). 

The variations of s{A) with A at different downstream locations x = i/2,i, and 2i are 
similar to each other: with A becoming more negative, i.e., further and further into the 
absolutely unstable regime the fluctuations sq and Sq become constant at levels that depend 
on the measuring location - the closer to the inlet where R becomes smaller the larger are 
the fluctuations. This behavior is reflecting the relation sq ~ oc that can be read off 
directly from Eq. ()3.5|) . 

The downstream reduction of the variance sq of the frequency fluctuations with increasing 
distance from the inlet and with increasing amplitude along the front is shown in Fig. [3 for 
the quintic GLE. There we compare the behavior of sq together with the front profiles of 
y^(P4py in the absolutely and in the convectively unstable regime close to the threshold 
A = for /i = 0.05. 

IV. CONCLUSION 

We have studied numerically the influence of small additive noise on pattern formation 
near a forwards and near an inverted bifurcation as described by a cubic and quintic GLE, 
respectively, when a finite group velocity v can blow the finite-amplitude part out of the 
system, i.e., in the vicinity of the so-called convective-absolute instability at A = v/vc-a{fJ') — 
1 = 0. The front that connects the inlet condition A{x = 0) = to the finite-amplitude 
downstream bulk part — R% is for the cubic GLE more sensitive to the applied noise 

strength than for the quintic case. This is partly related to the different magnitudes of the 
curvatures of the deterministic GLE potentials around the states A = and A = Rj^j: the 
resulting growth enhancement of fluctuations near A = is larger in the cubic than in the 
quintic case and in addition the damping of fluctuations near A = Rn is smaller in the cubic 
than in the quintic case. 

In the cubic case the transition between the regimes of small and large fluctuations of am- 
plitude, frequency, and wave number is shifted to a negative A into the absolutely unstable 
regime. Simultaneously the pattern growth length i{A) has there a characteristic inflection 
point that shows up as a peak in di/dA. In the quintic case all this occurs at the un- 
shifted convective-absolute threshold A = 0. Common to both cases is that the fluctuations 
decrease along the front in both regimes with growing pattern amplitude ■>/ < \A\'^ >. 
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For negative subcritical amplitude growth rates, < 0, we did not find noise-sustained, 
large-amplitude, backwards bifurcating patterns when A is positive: the nonlinear deter- 
ministic front dynamics of the quintic GLE blows any large-amplitude part downstream 
away from the inlet where A = and eventually any finite system is filled only with small- 
amplitude fluctuations of A around the stable fixed point A = of the unforced system. 
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APPENDIX 

Here we estimate the noise dependence of the downstream growth length i of the nonlinear 
structure in the convectively unstable regime of the cubic GLE where this structure is noise 
sustained. To that end we approximate i by the length where the mean squared amplitude 
Ciin{x) = {\Aiin{x)\'^) of the linear GLE has grown from the inlet value A{x = 0) = to, 
say, one half of the nonlinearly saturated bulk value (|v4p) ^ fi/2. So we solve the equation 

CiU^ = i) = ^fi (A.l) 

for i. Actually the linear solution may not hold there anymore. But as it will become 
obvious below the result is roughly independent of the coefficient chosen in Eq. ()A.1|) so also 
smaller numbers than ^ could be chosen here for a characteristic growth length. 

We evaluate the equal-time correlation Cunix) via the frequency integral of the spectrum 
Ciin{x,uj) of the time-displaced autocorrelation function of fluctuations of Aun at the same 
downstream position x. For large downstream distances x from the inlet this spectrum is 



given by 



with 

i^Q) = ±i\/lJ^c-a - fJ^-iuJ - iy/flc~a ■ (A. 3) 
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This spectrum f)A.2|) is strongly peaked at the center, c<j = 0, of the band of modes, 
—2y/fific-a < uj < 2^\x\ic-a-, that are amphfied in the convectively unstable regime. Thus, 
the aforementioned frequency integral may be approximated by 



Clin{x) 



7r-Clinix,Uj) ~ y/HlXc-aClin{x,UJ = 0) = CXp 



2iKi {lo=0)x 



(A.4) 



The last equality follows from Eq. ()A.2jl at = 0. Applying now the condition (jA.lj) one 
obtains 

1 

In 



21/2^3/4 



Using fic-a/ fJ' = (1 + in Eq. ()A.3|) one sees that iKi{uj = 0) = ^/JI 
for A -C 1 so that finally at fixed /x 

+ V2A + 0{A) \ (In a + const) . 



(A.5) 
(A.6) 
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FIG. 1: Convective-absolute instability boundaries (|2.1()|) for the unforced cubic and quintic GLE, 
respectively. For parameters below the respective curve front propagation is such that in the 
absence of noise the ^=0 state invades the A ^0 state. In the absolutely unstable parameter 
regime above the respective curve the ^=0 state recedes and the A 7^0 state expands (as long as 
the front is not hindered by a boundary). 
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R 

FIG. 2: Deformation of the front solution R{x) (a) of the deterministic quintic GLE by the 
boundary condition A{x = 0,t) = in the absolutely unstable regime for /x = 0.05 and v = 
(1 + A)vc-a as indicated. The spatial growth rate = dxlnRjx) (b) deviates from a freely 
propagating front with wave number g = that would show [2o| kl = | + \J~^^^J^ [thick dots in 
(b)] in the small-amplitude "linear" part of the front, (c) shows k versus R in comparison with the 
prediction ^\ ~ "71 ~ ^ stationary front in an infinite system at A = 0, /i = 0.05 

for which ii^=1.048. Thin dotted curves in (a)-(c) refer to a numerically obtained solution for 
A = at time 5 • 10^ which is not yet stationary. Here the profile is still moving to the right and 
in the absence of numerical "noise" we would expect this transient to approach the R = basic 
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FIG. 3: Reduced potentials V H2.13() corresponding to the real cubic (dashed line) and quintic 
(dotted line) GLE. For the cubic case the plot is independent of fJ-,g3- For the quintic case it 
depends on the combination jl = —gbfi/g'i- Here /2 = 0.05. 
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FIG. 4: Growth length £ (a) and its derivative d£/dA (b) versus A near the convective-absolute 

threshold for the cubic GLE at = 0.05 for various noise strengths a. Inset in (b) shows the 
variation of the peak of di/dA with a. 
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FIG. 5: Inverse of the standard deviations of frequency sq (a), wave number Sg (b), and amplitude 
modulus sr (c) for the stochastic cubic GLE. Results are reduced by S (|,'-{.()|) or ^y^~\A\^~>, 
respectively, and plotted as functions of A for three downstream locations x = i/2, i, and 2i. (d) 



shows the growth length £ of < \A\'^ > together with its derivative di/dA. Piecewise straight 
lines are guides to the eye. Parameters are /U = 0.05 and a = 2.5 • 10~^. 
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FIG. 6: Inverse of the standard deviations of frequency sq (a), wave number Sg (b), and amplitude 
modulus sr (c) for the stochastic quintic GLE. Results are reduced by S or a/ < >, 

respectively, and plotted as functions of A for three downstream locations x = i/2, i, and 2i. (d) 



shows the growth length £ of < \A\'^ > together with its derivative di/dA. Piecewise straight 
lines are guides to the eye. Parameters are /U = 0.05 and a = 10~^. 
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FIG. 7: Spatial variation of the standard deviation sq of the frequency reduced by S (a) and of 
a/ < l^p > (b) for the quintic GLE in the absolutely and convectively regime at A = —0.1273 and 
A = 0.0908, respectively. After the integration time of T = 5 • 10^ used in this plot sq was for 



A = 0.0908 not yet fully stationary in the growth region of ^/ < \A\'^ >. Parameters are /j, = 0.05 
and a = 10'^. 



